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1 .  Short  Description  of  the  Lidar  System  Used  for  this 
Research,  Improvements  and  Modifications 


The  lidar  used  for  this  project  is  a  mobile  system  with  a 
ruby  and  a  neodymium  glass  laser  transmitter.  Both  lasers 
are  operated  with  and  without  frequency  doubling,  resulting 
in  available  frequencies  of  347,  530,  694,  and  1060  nm.  The 
backscattered  light  is  collected  by  a  30  cm  dia.  reflecting 
telescope,  and  detected  by  an  EMI  9813  photomultiplier  tube 
with  bialkali  cathode  (347  and  530  nm)  and  a  large-area  PIN- 
diode  with  broadband  amplifier  (694  and  1060  nm) .  The  signals 
are  digitized  by  a  Biomation  8100  transient  recorder  with  8 
bit  amplitude  resolution,  10  ns  minimum  sample  interval  and 
2000  bytes  memory,  and  further  processed  by  an  on-line  com¬ 
puter.  The  range-corrected  and  energy-normalized  backscatter 
profiles  are  represented  on  a  display  tube  and  plotted  by  an 
X-Y  plotter,  both  incorporated  in  the  lidar  system.  The  dig¬ 
itized  original  returns  are  stored  on  floppy  disks  or  mag¬ 
netic  tape  cassette,  together  with  housekeeping  and  other 
auxiliary  data.  ... 


In  the  past,  the  change  from  the  fundamental  wavelength  to 
the  first  harmonic  was  done  for  both  lasers  through  manual 
replacement  of  the  beam  expanding  telescopes  which  are  focused 
and  aligned  respectively  for  only  one  wavelength.  The  fre¬ 
quency  doublers  were  connected  to  one  of  the  telescopes  each. 
In  order  to  avoid  this  time-consuming  procedure,  at  least  for 
the  neodymium  laser,  an  achromatic  telescope  with  common  focus 
for  530  and  1060  has  been  calculated  and  constructed,  which 
can  remain  permanently  in  the  beam  path,  so  that  the  neodymium 


laser  transmits  both  wavelengths  simultaneously.  The  time  lag 
between  measurement  in  both  wavelengths  is  thus  only  limited 
by  the  data  processing  time,  which  is  a  few  seconds  if  a  real¬ 
time  display  and  plot  of  the  normalized  profiles  is  dispensed 
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Unfortunately,  a  similar  procedure  is  not  applicable  to  the 
ruby  laser  because  no  UV-permeable  optical  glasses  with  suf¬ 
ficiently  different  dispersion  are  available.  For  the  present 
task  the  347  nm  wavelength  is  less  important  anyway.  Never¬ 
theless,  an  optical  construction  is  in  progress  allowing  the 
expanding  of  the  694  nm  beam  by  the  new  achromatic  telescope, 
too.  For  this  purpose,  a  beam  divider  and  a  deflection  mirror 
will  be  built  into  the  system. 

The  originally  fixed  gain  of  the  diode  detector  amplifier 
proved  to  be  too  high  for  the  recording  of  the  strong  back- 
scatter  signals  from  clouds  and  fog  without  overload  effects. 
For  attenuation  of  those  signals  two  neutral  filters  with 
transmissivities  of  50  and  25%,  respectively,  were  inserted 
into  free  positions  of  the  interference  wheel.  Furthermore, 
by  switching  to  a  smaller  feedback  resistor,  it  is  now  pos¬ 
sible  to  reduce  the  gain  by  a  factor  of  11.6. 


Program  Developments 


The  application  of  Klett's  evaluation  method  for  lidar  returns 
requires  a  considerable  amount  of  computation  work.  To  do  this, 
some  time  ago  a  program  was  developed  for  the  lidar  on-line 
computer.  Profiles  of  the  extinction  coefficient  calculated 
from  the  lidar  data  by  Klett’s  method  and  the  "standard"  method, 
where  the  boundary  value  is  chosen  close  to  the  lidar  system 
and  the  integration  of  the  lidar  equation  proceeds  away  from 
it  (called  "CRBV  method"  in  the  previous  interim  report) , 
were  plotted  by  the  above-mentioned  X-Y-plotter.  This  plotter, 
however,  is  not  very  comfortable  and  precise,  especially  if 
curves  with  heavy  fluctuations  are  to  be  displayed.  Further¬ 
more,  the  memory  capacity  of  the  lidar  on-line  computer 
(32  k  bytes)  is  somewhat  too  small  for  the  handling  of  large 
amounts  of  data  like  digitized  lidar  returns.  For  these  reasons 
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we  decided  to  transfer  all  raw  lidar  data  via  magnetic  tape 
to  the  larger  computer  which  is  available  since  about  a  year. 
This  computer  offers  512  k  RA.'  and  400  megabytes  disc  memory. 
The  operating  system  can  support  up  to  63  simultaneous  users. 
Each  user  is  totally  independent  and  has  768  kilobytes  vir¬ 
tual  memory  at  his  disposal.  An  extensive  library  of  auxili¬ 
ary  programs  and  subroutines  is  available,  for  instance, 
routines  for  the  precise  display  of  curves  on  a  graphic 
terminal  screen  or  by  an  incremental  plotter. 

The  programs  for  the  application  of  the  Klett  and  standard 
method  for  the  evaluation  of  lidar  returns  have  been  rewrit¬ 
ten  in  ECKERAN  language  for  this  new  computer.  Furthermore, 
extensions  have  been  made  for  performing  the  following  tasks: 

-  listing  of  the  digitized  original  lidar  return  as  function 
of  distance 

-  linear  plot  of  the  original,  not  range-corrected  lidar 
return 

-  linear  plot  of  the  range-corrected  lidar  profile  in 
absolute  units 

-  logarithmic  plot  of  the  range-corrected  lidar  profile, 

-  linear  plot  of  the  extinction  coefficients  and  optical 
depth  profiles,  calculated  by  the  standard  and  Klett’ s 
method 

-  listing  of  the  extinction  coefficient  and  optical  depth 
versus  distance 

-  storage  of  digitized  raw  lidar  data  and  calculated 
extinction  coefficient  and  optical  depth  profiles  on 
magnetic  tape  in  any  desired  format. 

The  disc  memory  capacity  of  the  computer  is  large  enough  for 
permanent  storage  of  numerous  lidar  returns  and  calculated 
profiles  of  optical  parameters,  which  greatly  facilitates 
the  computation  work.  Examples  of  listed  and  plotted  profiles 


will  be  described  in  Chapter  4. 
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3.  Measurements 


In  general,  situations  with  strong  and  inhomogeneous  fog  or 
aerosol  backscatter  in  the  near-ground  atmosphere  occur  less 
frequently  during  summertime.  However,  in  some  cases  we  were 
able  to  get  lidar  returns  from  a  turbid  atmosphere,  with  scat¬ 
tering  coefficients  distributed  both  homogeneously  and  inhomo- 
geneously  along  the  laser  beam  path  (measurements  before  April 
1983  have  already  been  listed  in  the  previous  interim  reports). 

On  June  29,  a  very  humid  air  mass  with  high  aerosol  concentra¬ 
tion  was  present,  and  in  the  morning  some  low-level  cumulus 
fractus  clouds  occurred.  Spatial  inhomogeneities  of  the  back¬ 
scatter  coefficient  were  observed  also  outside  the  cloud, 
caused  by  partial  precondensation  (growth  of  the  aerosol  par¬ 
ticles  due  to  water  absorption).  Towards  noon,  the  humidity 
decreased,  and  the  backscatter  became  homogeneous  in  the  lowest 
1000  m  of  the  atmosphere.  Five  measurement  series  were  obtained 
in  the  wavelengths  530,  694,  and  1060  nm. 

0n  July  14,  similar  aerosol  conditions  occurred  after  a  thunder 
shower.  Again,  cumulus  fractus  clouds  were  observed  about  400  m 
above  ground  level  during  the  morning  hours,  and  strong  and  ini 
tially  inhomogeneous  aerosol  backscatter  below.  Like  on  June  29 
the  backscatter  homogenized  with  decreasing  humidity  towards 
noon  time,  but  the  concentration  of  the  dry  aerosol  was  consid¬ 
erably  higher  in  this  case.  Eleven  series  of  measurements  were 
performed  in  each  wavelength.  One  of  them  will  be  discussed  in 
more  detail  in  the  next  section  of  this  report. 

During  the  next  day,  July  1 5 ,  the  aerosol  constitution  remained 
nearly  unchanged.  In  the  afternoon  some  measurements  were  per¬ 
formed  under  homogeneous  conditions. 

On  August  1 2 ,  a  shallow  layer  of  low-concentrated  ground  fog 
was  present  during  the  late  night  and  early  morning  hours. 

This  layer,  however,  dissolved  then  very  rapidly,  and  so  only 
one  lidar  return  was  obtained  at  6  AM  CET. 


4.  Further  Example  of  Extinction  Coefficient  Profiles 
Calculated  by  the  Klett  and  CRBV  Method,  Application 
of  Aerosol  Data  for  Derivation  of  Klett  Boundary  Values 

The  basic  principle  of  Klett' s  method  is  to  choose  the  boundary 
value  for  the  extinction  coefficient  at  the  far  end  of  the  range 
interval,  not  at  the  close  one.  The  dependence  of  this  boundary 
value  (BV)  upon  the  resulting  extinction  coefficient  profile 
is  the  less  the  higher  the  optical  depth  is  at  that  point.  In 
the  previous  interim  report  from  21  March  we  pointed  out  that 
under  dry  aerosol  conditions ,  even  at  extremely  high  concentra¬ 
tions,  the  optical  depth  at  a  3  km  range  is  in  most  cases  not 
high  enough  to  keep  the  influence  of  the  BV  sufficiently  low. 
Increasing  the  lidar  range  is  not  always  applicable  due  to  the 
signal-to-noise  ratio,  rapidly  increasing  with  distance.  In 
such  cases,  therefore,  a  BV  more  or  less  close  to  reality 
must  be  found.  If  the  spatial  aerosol  distribution  is  suffi¬ 
ciently  homogeneous,  the  slope  or  average  slope  of  the  range- 
corrected  backscatter  profile  may  be  used  for  this  purpose. 

In  our  case,  we  have  the  advantage  to  operate  an  unmanned 
measuring  station  for  meteorological  and  aerosol  data  at  a 
mountain  peak  at  1650  m  altitude  above  MSL,  2.55  km  apart  from 
the  Institute,  in  connection  with  an  atmospheric  transmission 
path.  The  station  can  be  switched  on  and  off  by  a  radio  signal. 
Aerosol  size  distributions  are  measured  by  a  Knollenberg  clas¬ 
sical  scattering  aerosol  spectrometer.  Model  PDS-200.  From  the 
measured  particle  size  spectra,  extinction  coefficients  as  a 
function  of  wavelength  can  be  calculated  using  van  de  Hulst's 
approximation  formula.  On  July  14,  the  particle  spectrometer 
was  operated  during  the  lidar  measurements,  and  extinction 
coefficients  were  calculated  for  the  lidar  wavelengths. 

Fig.  1  presents  an  original  lidar  return  in  694  nm  (ruby  laser) 
wavelength,  normalized  fcith  respect  to  transient  recorder  in¬ 
put  range  and  laser  energy  but  not  to  range ,  i . e .  the  curve  is 
proportional  to  the  electrical  signal  from  the  photodetector 


6 


(the  peak  near  the  end  of  the  range  is  the  reflex  from  a  moun¬ 
tain  face) .  The  profile  was  calculated  by  the  new  program 
mentioned  in  Chapter  2,  displayed  on  a  graphic  terminal  screen 
and  hardcopied.  At  the  upper  left  corner  are  given  the  profile 
identifier,  composed  of  a  two-digit  series  number,  a  wavelength 
code  (A) ,  and  the  date,  and,  below  it,  the  maximum  abscissa 
(range)  and  ordinate  values.  The  ordinate  units  here  are 
volts/Joule.  A  listing  of  the  digitized  profile  amplitude 
as  function  of  range  in  meters  is  presented  in  Table  1  in  the 
appendix.  In  this  table  the  data  proceed  linewise  (horizontally) 
with  distance.  Each  line  contains  six  pairs  of  distance  (in  m) 
and  amplitude  values;  the  amplitudes  in  this  case  are  the  dig¬ 
ital  numbers  delivered  by  the  transient  recorder  (or  the  sum 
of  them,  if  more  than  one  signal  is  averaged)  minus  the  back¬ 
ground  level,  without  any  normalization  factor.  The  format  of 
the  lines  is  described  by  the  FORTRAN- statement  FORMAT (6(F7.1 ,15) ) . 
In  this  presentation  mode  the  limited  amplitude  resolution  of 
the  signal  digitizer  remains  clearly  evident.  The  whole  table 
is  stored  as  a  data  file  in  the  disc  memory  of  the  computer 
in  the  same  format,  i.e.  one  line  with  six  data  pairs  is  one 
record.  The  transfer  to  other  carriers,  e.g.  magnetic  tape,  is 
thus  very  simple  by  calling  a  special  service  program.  The 
example  is  a  Ildar  return  from  a  homogeneous  atmosphere;  a 
return  from  an  inhomogeneous  cloud  is  listed  in  Table  2,  it 
is  the  example  discussed  in  the  1st  interim  report  from  Feb¬ 
ruary  10,  1983. 

Fig.  2  shows  the  range  corrected  backscatter  profile  calculated 
from  the  lidar  return  Fig.  1;  since  calibration  factors  are 
known  for  each  wavelength,  the  profile  amplitude,  presenting 
the  product  B*t2  ( 8*backscatter  coefficient,  x*transmissivity 
factor),  can  be  given  in  absolute  units,  m-1  sterad""1 .  In  the 
upper  left  corner  again  are  shown  the  identifier  and  the  ab- 


scissa  and  ordinate  maxima,  XM  and  YM,  respectively. 
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For  some  purposes  a  logarithmic  plot  of  the  range-corrected 
profile  is  desirable.  In  Fig.  3  our  profile  is  presented  in 
such  a  wav,  where  the  decadic  logarithm  has  been  chosen.  The 
full  ordinate  range  is  one  decade  in  linear  scale  or  unity 
in  logarithmic  scale,  as  indicated  in  the  upper  left  corner. 

The  nearly  linear  decay  of  the  lidar  profile  in  this  presen¬ 
tation  indicated  almost  homogeneous  aerosol  conditions,  and 
an  average  extinction  coefficient  .315  km-1  is  easily  derived. 

As  already  mentioned  above,  for  our  example  we  have  aerosol 
concentration  and  si2e  distribution  data  from  the  far  end  of 
the  lidar  range  interval,  which  can  be  used  for  the  derivation 
of  the  Klett  boundary  value.  Fig.  4  shows  hourly  averaged  par¬ 
ticle  size  distribution  curves  measured  near  the  end  of  the 
lidar  range  by  an  optical  sprectrometer .  The  size  range  is 
.32  to  20  microns  in  diameter.  From  the  spectrum  closest  to 
the  time  of  the  lidar  measurement,  van  de  Hulst's  approxima¬ 
tion  results  in  an  extinction  coefficient  of  1.57  km-1  for 
the  ruby  wavelength,  where  a  refractive  index  of  1.419-.0331i 
was  assumed ,  according  to  Shettle  a  Fenn ’ s  urban  aerosol  model 
at  a  relative  humidity  of  80%  as  measured.  If  we  use  the  ex¬ 
tinction  coefficient  just  mentioned  and  apply  Klett's  method 
to  our  lidar  profile,  we  get  the  extinction  coefficient  profile 
shown  by  the  solid  curve  in  Fig.  5.  The  dashed  line  is  the 
optical  depth  as  function  of  range.  Above  left  are  given  the 
file  identifier  as  before,  the  maximum  abscissa  value  (3000  m 
in  this  case) ,  the  ordinate  maxima  for  extinction  coefficient 
(YMS)  and  the  optical  depth  (YMT) ,  the  chosen  value  of  the 
exponent  Q  in  the  exponential  relation  between  backscatter 
and  extinction  coefficients  (here:  1),  the  reference  distance 
REFD ,  and  the  boundary  value  SR. 

We  see  in  Fig.  5  that  the  BV  calculated  from  the  measured 
aerosol  data  is  obviously  too  high,  as  indicated  by  the  unrea¬ 
listic  increase  of  the  extinction  coefficient  towards  the  end 


of  the  interval.  Actually  the  sprectrometer-derived  extinction 
is  almost  five  times  higher  than  the  slope  value.  The  reason 
for  this  discrepancy,  which  was  observed  in  other  cases,  too, 
is  not  yet  known.  A  local  enhancement  of  particle  size  at  the 
spectrometer  site  due  to  cloud  formation  or  precondensation 
is  improbable,  since  the  measured  relative  humidity  never  ex¬ 
ceeded  80%.  Nevertheless,  the  extinction  coefficient  derived 
from  aerosol  data  assessed  at  ground  level  is  much  smaller 
than  that  from  the  rear  end  of  the  lidar  range  interval, 
namely  .185  km  1 ,  calculated  from  optical  spectrometer  data 
and  .175  km  1 ,  calculated  from  impactor  data.  These  values 
are  even  smaller  than  the  extinction  coefficient  derived  from 
the  lidar  profile  slope.  Further  investigations  seem  to  be 
required  in  the  future  to  find  the  reason  for  these  discrep¬ 
ancies  . 

If  the  Klett  method  is  applied  with  the  slope  extinction  co¬ 
efficient  as  BV,  the  profile  presented  in  Fig.  6  is  obtained. 

As  expected,  a  nearly  constant  extinction  coefficient  results. 
The  same  is  true  if  the  standard  or  CRBV  method  is  applied 
with  that  boundary  value  at  the  near  end  of  the  interval,  as 
shown  in  Fig.  7.  The  meaning  of  the  solid  and  dashed  curves 
is  the  same  as  before.  Above  left  are  printed  the  identifier, 
the  close  and  far  ends  of  the  range  interval,  XR  and  XE,  the 
coordinate  maxima  for  range,  extinction  coefficient  and  opti¬ 
cal  depth,  XM,  YMS  and  YNT,  respectively,  the  exponent  Q,  and 
the  boundaryvalue  SR,  here  valid  for  range  XR. 

A  comparison  of  Figs.  6  and  7  suggests  both  methods  to  be  equi¬ 
valent  for  lidar  returns  from  a  homogeneous  atmosphere.  However, 
the  Klett  method  is  still  superior  in  so  far  that  the  choice 
of  the  BV  is  much  less  critical. 

The  extinction  and  optical  depth  profiles  obtained  by  the  Klett 
and  CRBV  methods  may  be  stored  permanently  as  data  files  in  the 
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disc  memory  of  the  computer.  Since  these  files  are  very  ex¬ 
tensive,  we  present  in  the  appendix.  Table  3,  only  the  begin¬ 
ning  of  such  a  file,  here,  the  Klett  profile  shown  in  Fig.  5. 

The  first  line  contains  the  maximum  number  of  profile  data 
points  and  the  maximum  extinction  coefficient  and  optical  depth, 
the  second  one  the  reference  distance  in  m,  the  boundary  value 
in  m-1 ,  and  the  exponent  Q.  The  lines  following  contain  the 
range  in  m,  the  extinction  coefficient  in  m  1 ,  and  the  optical 
depth . 


5 .  Future  Plans 

The  soft-ware  is  now  complete  except  some  improvements  like 
plotting  of  coordinate  scales  in  the  diagrams.  One  of  the  tasks 
of  this  contract  work  is  the  examination  of  the  appropriateness 
of  the  relationship 

fl(r)  =  a*o(r)q  . 

Only  if  this  relationship  is  valid  an  analytical  solution  of 
the  lidar  equation  exists.  The  examination  of  the  validity  is 
somewhat  difficult.  A  possible  way  to  do  this  is  to  correlate 
the  absolute  lidar  amplitudes  close  to  the  lidar,  where  the 
extinction  can  be  neglected,  with  extinction  coefficients  at 
the  same  range,  calculated  by  any  appropriate  method.  For  this 
purpose  as  many  as  possible  measurements  are  required,  during 
fog  and  cloud  conditions  from  within  the  fog.  This  is  neces¬ 
sary  to  get  absolute  fog  backscatter  data .  In  order  to  get 
such  data  it  will  be  necessary  to  go  with  the  lidar  to  sites 
where  ground  fog  is  more  frequent  than  at  our  Institute.  Such 
places  exist  in  the  plains  outside  of  our  valley;  e.g.  in  the 
"Murnauer  Moor"  or  in  other  areas  not  too  far  from  here. 

Since  only  a  part  of  such  excursions  can  be  expected  to  be 
fully  successful ,  it  is  doubtful  that  this  task  on  the  basis 
of  experiments,  data  evaluation,  and  theoretical  work  can  be 
finished  during  the  current  contract. 


On  the  other  hand,  the  lidar  shots  close  to  the  aerosol 
spectrometer  at  the  mountain  station  should  be  continued 
to  find  out  whether  a  general  discrepancy  exists  between 
spectrometer-  and  lidar-derived  extinction  coefficients. 


Garmisch-Partenkirchen ,  September  8,  1983 


(Dr.  R.  Reiter) 
Director 
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